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ABSTRACT 

It is assumed that accreting neutron stars in low-mass X-ray binaries are heated due 
'. to the compression of the existing crust by the freshly accreted matter which gives 

O j rise to a variety of nuclear reactions in the crust. It has been shown that most of 

the energy is released deep in the crust by pycnonuclear reactions involving low-Z 
elements (the deep-crustal heating scenario). In this paper we discuss if neutron stars 
in so-called very-faint X-ray transients (VFXTs; those which have outburst peak 2- 
10 keV X-ray luminosities < 1 x 10 36 erg s _1 ) can be used to test this deep-crustal 
heating model. Unfortunately we cannot conclusively answer this question because of 
the large uncertainties in our estimates of the accretion rate history of those VFXTs, 
both the short-term (less than a few tens of thousands of years) and the one throughout 
their lifetime. The latter is particularly important because it can be so low that the 
neutron stars might not have accreted enough matter to become massive enough that 
enhanced core cooling processes become active. Therefore, they could be relatively 
warm compared to other systems for which such enhanced cooling processed have 
been inferred. However, the amount of matter can also not be too low because then 
the crust might not have been replaced significantly by accreted matter and thus a 
hybrid crust of partly accreted and partly original, albeit further compressed matter, 
might be present. This would inhibit the full range of pycnonuclear reactions to occur 
and therefore very likely decreasing the amount of heat deposited in the crust. More 
detailed calculations of the heating and cooling properties of such hybrid crusts have 
to be performed to be conclusive. Furthermore, better understanding is needed how a 
hybrid crust affects other properties such as the thermal conductivity. We also show 
that some individual neutron-star low-mass X-ray binaries might have hybrid crusts as 
well as possibly many of the neutron stars in high-mass X-ray binaries. This has to be 
taken into account when studying the cooling properties of those systems when they 
are in quiescence. In addition, we show that the VFXTs are likely not the dominate 
transients that are associated with the brightest (~ 10 33 erg s _1 ) low-luminosity X-ray 
sources in globular clusters as was previously hypothesized. 
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1 INTRODUCTION 

Neutron-star low mass X-ray binaries (LMXBs) harbor neu- 
tron stars which are accreting matter from a close-by low- 
mass (typically < 1 Mq) companion star which transfer 
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mass due to Roche-lobe overflow. In most systems the neu- 
tron star does not continuously accrete matter. Generally, 
those systems (called X-ray transients) are in their quiescent 
state in which they do not accrete at all or only at a very low 
rate. Only occasionally they exhibit bright X-ray outbursts 
during which their X-ray luminosities increase by several or- 
ders of magnitude. Such outbursts are most probably caused 
by a very large increase in the mass accretion rate onto the 
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neutron stars due to insta bilities in the accretion disc (see 
the review bv lLasotallioOlb . The X-ray transients can be di- 
vided in sub-groups based on their peak 2-10 keV X-ray lu- 
minosities in outburst. In particular, a special sub-group has 
been called the very-faint X-ray transients (VFXTs) that 
have pe ak luminosities betwee n 1 x 10 34 and 1 x 10 36 erg 
s _1 (see IWiinands et al.l 120061 ) . whereas the brighter tran- 
sients have peak luminosities of io 36-39 erg s _1 . 

In quiescence, neutron-star transients can still be de- 
tected using sensitive X-ray satellites, and it has been 
found that in many systems a soft, most likely thermal, 
component is present with a ty pical black-body temper 



ature of 0.1-0.3 keV (see, e.g. , van Paradiis et al., 
lAsai et al. 1996; Ca mpana et al.ll 19981 ; iRutledge etal\ 



1987 



1999 



and references to those papers). In addition, for many 
systems an additional spectral component above 2 keV 
has also been det ected (the non-t h ermal power-law com - 
ponent; see, e.g., lAsai et all 1 19961 ; IRutledge et~aH Il999l) . 
which can even dominate the 0.5-10 keV X-ray flux in some 
systems (e.g., Campana et al.l l2002al; Ijonker et al.1 12004 



IWiinands et al.ll2005l ; iDegenaar et alfe^bf l The origin of 
the non-th ermal component is not clear (see, e.g., the dis- 
cussions in lCampana et a l. 1998; D egenaar et alj|2012bf ) but 
it is generally assumed that the soft component is the ther- 
mal emission from the neutron star surface, either due to 
very-low level residual accretion onto the surface or due to 
the cooling of the neutron star that has been heated by the 
matter accreted during the outbursts. 

During the accretion phases, matter accumulates on the 
surface of the neutron star. This matter compresses the un- 
derlying layers of the neutron star crust. If the accretion 
continues long enough the original catalyzed crust can be 
co mpletely replaced by a new crust ma de of accreted mat- 
ter (|Satdll979l lHaensel fc ZdunikHl990bl ). The original crust 
is pushed down into the neutron star until it fuses together 
with the core. The composition of the accreted crust should 
be quite different from t he original, catalyzed cru st, i.e., 
richer in low-Z elements (|Haensel fc Zdunikl 1 1990bh . It has 
been postulated that when the accreted matter sinks into the 
crust due to the compression induced by freshly accreted ma- 
terial onto the star, a chain of non-equilibrium reactions oc- 
curs in the crust that generates heat (elec tron captures, neu- 
tron drips and pycnonuclear reaction s; lHaensel fc Zdunikl 
Il990al . 120031 . 120081 ; iGupta et alj|2007h . Most of the heat is 
released deep in the crust (at densities > 10 12 g cm -1 ) due 
to pycnonuclear reactions involving low-Z elements. This 
heat is conducted inwards, heating the core, and outwards, 
where it is emitted as thermal emission from the surface. 
This model has bee n called the "deep crustal heating model" 
(|Brown et al.ll 19981 ). 

In Section [5] we briefly describe the deep crustal heat- 
ing model and compare the model with the available data. 
We also calculate the time-scale on which the core reacts 
to changes in the long-term averaged accretion rate. In 
Section [3] we calculate, in the frame work of the deep- 
crustal heating model, the expected quiescent luminosity of 
neutron-star VFXTs, in order to use those systems to test 
the deep-crustal heating model. We argue that it might be 
possible that during their life those systems might not have 
accreted enough matter to have fully replaced their orig- 
inal neutron-star crust with an accreted one which could 
significantly inhibit the pycnonuclear heating reactions. In 



Section [4] we discuss how the VFXTs can still be used to 
test the model and also discuss potential other sources which 
might harbor neutron stars with only partly accreted crusts. 



2 QUIESCENT NEUTRON STARS AND THE 
DEEP CRUSTAL HEATING MODEL 

In this heating/cooling model one can calculate the ther- 
mal state of the neutron star with a simple energy balance 
consideration by writing 



dE th 



dT 

C'v —r- = H — L-, 
dt dt 1 



(1) 



where Eth is the thermal energy of the neutron star, Gy 
its total specific heat, and T its core temperature. H is the 
total heating rate, and the two energy sinks are the star's 
thermal photon luminosity, L 7 , and its neutrino luminosity 
L v . 

In the deep crustal heating scenario H is taken as a time 
average 



H 



10" 



(M) 



<M>- 



10" 11 M Q yr- 1 1.5 MeV 



nuc — 1 

erg s 



(2) 



where (M) is the long term time-averaged mass accretion 
rate onto the neutron star, Q n uc the amount of heat, per 
accreted nucleon, deposited in the crust and m u the atomic 
mass unit. Theoretical predictions (e.g., lHaensel fc Zdunikl 
l2008f ) obtain values for Q nU c between 1 and 2 MeV. 

When (M) has been stable for a long enough time the 
neutron star is in thermal equilibrium (see Section 12.21 for 
estimates of the thermal response time-scale). Hence, from 
Eq. {TJ with dT /dt = 0, one obtains the expected L 7 , or the 
observable quiescent thermal luminosity iJ3 > as 

L q = (H) - (U) . (3) 

I Brown et a l. (1998) showed that when (L„) is negligible the 
very simple relation 



L q = (H) = (M)Q nuc /m u 



(4) 



agreed with several observa tions of qui e scent neutron-star 
LMXB transients. However, IColpi et al.l (|200ll ) emphasized 
that many such systems are hot enough that (L„) is not 
necessarily negligible and could even be in a regime where 
L v 3> L-y . In the case that fast neutrino emission is possible 
(as one would expect for a massive neutron star) this would 
explain the low values of L q observed from several systems 
that are discrepant with the simplified Eq. Q since with a 
very large L v one could have L q <C (M)Q nuc /m u . 

Neutrino emission processes can be roughly divided in 
two categories, either "slow" or "fast" that differ in their 
temp e rature dependence and efficiency (|Y akovlev fc Pethickl 



12004 iPage et al.1 120061 ; iPage fc Reddvl 120061 ). Slow pro- 
cesses include the modified Urea ("MUrca") and sev- 
eral bremsstrahlung processes whose emissivities can be 



roughly written as e s v 



QTg erg cm -3 s" 1 



where 



1 We will designate the photon luminosity by L 7 or L q , using 
the former when considering it as an energy sink, as in Eq. ((TJ, 
and the latter when it is seen as an observed thermal luminosity. 
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T 9 = T/10 9 K. Values for Q range from ~ 10 19 for 
the bremsstrahlung processes to 10 21 for the MUrca. A 
higher emissivity for the MUrca has been proposed by 
IVoskrensenskv fc Senatorovl (|l986l ), called the Medium- 
Modified Urea process ("MMUrca"), resulting from the soft- 
ening of the pion mode and that is a precursor to the pion 
condensate. The MMUrca has a Q that reach 10 23 and it 
smoothly merges into the pion condensate emissivity (see 
below) with growing density. The fast processes have emis- 
sivities 4f st = QTcj erg cm -3 s _1 . The simplest and most 
efficient fast process is the direct Urea ("DUrca") with nu- 
cleons and this process has Q ~ 10 27 . In the presence of 
hyperons other DUrca processes are possible with slightly 
reduced efficiencies. Other fast, but less efficient, processes 
are possible in the presence of a meson condensate (either 
pion or kaon) that also have a T 6 dependence and Q in the 
range 10 24 - 10 26 for pions and 10 23 - 10 25 for kaons. Finally, 
if the neutron star has an inner core with deconfined quark 
matter, similar DUrca processes are present with efficiencies 
that can match the nucleon DUrca one. 

The resulting neutrino luminosity is then given by 

LsloW 4 3 ^slow^ ArsloW rjiS /r\ 

~ -ttR Q T g =N T 9 (5) 

and 

L* ast = ^ttR 3 p ■ Q fast T 9 6 = iV fast T 9 6 (6) 

where R is the radius of the neutron star and R p the radius 
of the inner core (i.e., high density) region where the given 
fast process is acting. There are theoretical uncertainties on 
Q of a factor of a few for the slow processes and the DUrca 
ones, but they are much larger in the case of the meson 
condensates. Moreover, for the slow processes R ~ 10 km 
while for the fast ones R p can range from ~ km, i.e. almost 
no fast neutrino emission, to almost 10 km depending on the 
equation of state (in the case it allows such processes) and 
the mass of the neutron star. We take as a typical value 
R p ~ 5 km but varying it can change JV fast by almost three 
orders of magnitude. 

In the presence of pairing (causing superfluidity and/or 
superconductivity), neutri no emission pro cesses can be 
strongly altered (see, e.g.. iPage et alj|2006l ). In low mass 
neutron stars where the MUrca processes can be strongly 
suppressed by pairing, bremsstrahlung processes may be- 
come the dominant sources of neutrinos. In the fast cooling 
scenarios L^ st can also be significantly reduced by pairing. 
Moreover, pairing opens a new neutrino emission channel 
from the constant formation and breaking of Cooper pairs 
(dubbed the "pair breaking and formation" or "PBF" pro- 
cess). The corresponding emissivity can be roughly approxi- 
mated by e PBF m QTg erg cm -3 s _1 where Q can reach 10 22 
in optimal conditions depending on the t ype of pairing an d 
its corresponding critical temperature T c (|Page et al.ll2009l ). 

The photon luminosity is simply L 1 = AirR^asBiT^ ) 4 
where <jsb is the Stefan-Boltzmann constant and T e °° the 
star's red-shifted effective temperature. T e has to be related 
to the internal temperature T and, as a simple rule, one can 
use T e « 10 6 T g 1/2 K. This implies that L 7 ~ 7 x 10 32 T g 2 
erg s~ . A detailed study of accreted neutron star en- 
velopes and the result ing T e — T relationship can b e found in 
IPotekhin et~aH (|l99it ) and lYakovlev et all \2004 ). The lat- 
ter work shows that the T dependence of L 1 range between 

© 0000 RAS, MNRAS OOO.rTlfMl 



T 1 ' 7 to T 2 ' 3 , depending on the actual chemical composition 
of the accreted envelope. 

2.1 Comparison of deep crustal heating with data 

Predictions for L q as a function of (M) for the various neu- 
trino emission scenarios described above are shown in Fig. [T] 
and compared with data ( similar to what ha s been done b y 
lYakovlev fc Pethicklliobl ; iHeinke et ai]|2007l . l2009al . l20ld ) . 

In this figure the band "Heating" shows the average 
heating rate (H) from Eq. ((2j with Q n uc ranging from 1 to 2 
MeV. Any star located on this line balances its heating only 
by its L q and is thus in the photon cooling regime of Eq. @. 
However, an object located below it has significant neutrino 
losses: it is in the neutrino cooling regime and the difference 
between its observed L q and the corresponding value of (H) , 
at the same (M), on the "Heating" lines directly gives its 
L v from Eq. (©. 

Each curve in Fig. [T] for the various neutrino cool- 
ing scenarios, is given by the energy balance of Eq. ©. 
The para meters set we u s e is c onsistent with the one pro- 
posed by lYakovlev et ail (|2003l ). The photon luminosity is 
choosen as L 1 = 7 x 10 34 Tg erg s _1 and the heating rate, 
or luminosity, is obtained from Eq. Q with Q nllc = 1.5 
MeV. The four upper lines, "Brems.", "MUrca", "PBF", 
and "MMUrca" correspond to the four slow neutrino cool- 
ing scenarios. We use, for the correspondings L v in Eq. ©, 
JV Brems = 5 x 10 37 , jV MUl ' ca = 5 x 10 39 , iV PBF = 5 x 10 40 , 
and jY MMUrca = 5 x 10 41 . As mentioned above, the "PBF" 
scenario can have a very wide range of efficiencies and 
the case considered here corresponds to the most efficient 
one as deduced for maximal compatibility of the "minimal 
coolin g" scenario with data from isolated cooling neutron 
stars (|Page et al.l |2004| . 120091 ) , and with the interpretation 
of the observed cooling of the neutron sta r in the sup ernova 
remnant Cassiopeia A (|Page et al.ll201ll ; iPagd l2012f ). The 
"MMUrca" scenario actually covers a wide range of neu- 
trino emission efficiencies that strongly depend on the neu- 
tron star mass and smoothly merge into the pion cooling 
scenario. 

The three pairs of lines "Kaon" , "Pion" , and "DUrca" 
show the prediction for each corresponding scenario when 
maximum and strongly reduced L f * Bt are assumed. We use, 
for the correspondings L v in Eq. ((6]) , N Ksion ranging from 5 x 
10 40 to 5 x 10 43 , iV pion from 5 x 10 41 to 5 x 10 44 , and jV Durca 
from 5x 10 41 to 5x 10 45 . When L„ > L 1 one has Lf ast = (H) 
and since Ljf st cx T 6 , if the neutrino efficiency is reduced 
by a factor 10 3 , T 2 must be 10 times higher for Ljf st to 
keep matching (H). Since L 7 oc T 2 , the resulting predicted 
L q is an order of magnitude higher. Reduced efficiency of 
Ljf st can be either due to an emissivity lower than quoted 
above or due to a smaller R p . Thus, to each one of these fast 
neutrino processes corresponds a band of at least one order 
of magnitude width in predicted L q . 

Comparison of these predictions with the data first 
shows that only a few objects are on the "photon cooling" 
line, i.e., are described by Eq. @, and that most observed 
quiescent neutron-star LMXBs have an in ferred L g which re - 
quires significant neutrino emission (e.g.. lColpi et al]|200ll ). 
In most case one sees that L v dominates over L 1 by one to 
two orders of magnitude. Enhanced neutrino emission im- 
plies a large neutron star mass that must have been accreted 
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(M) (M@ yr- 1 ) 

Figure 1. Locations of steady state quiescent luminosity L q ver- 
sus averaged mass acretion rate, (M) for a series of cooling scenar- 
ios. Each curve, for the various labeled neutrino cooling scenar- 
ios, is given by the energy balance of Eq. (J3j and the "Heating" 
band show the predicted range of (H) values: (see Section 12.11 
for details). Dotted lines are tracks of constant Tthi from Eq. (SJ, 
for high specific heat and, in parenthesis, low specific heat (i.e, 
CV = 10 30 T erg K" 1 and 10 29 T erg K -1 , or f SF = 1 and 0.1, 
respectively). The displayed observational data are taken from 
iHeinke et at] \2Q1(\) . 



during the lifetime of the X-ray binary, unless most neutron 
stars can be born this massive. The latter possibility would 
conflict with the observation of young isolated cooling neu- 
tron stars: the cooling of these stars is driven by neutrino 
emission and their observed luminosity can be reproduced 
with a L v simila r to the "PBF" model of Fig. [1] (see, e.g., 
iPage et al.ll2009l fl. 

The conclusion is that, unless neutron stars in LMXB 
are born more massive that isolated ones, most neutron stars 
shown in Fig. [T] must have accreted a significant amount of 
mass during their lifetime. A way around this conclusion is 
that the presently deduced (M) is not representative of the 
recent history of these systems (i.e., (M) was significantly 
lower in the past) and for this we examine in the next sub- 
section their thermal inertia. An alternative is that there is 
something missing in the deep crustal heating model. 



2.2 Thermal response time-scale 

Since human X-ray observations of neutron-star LMXBs 
only span, in the best case, a few decades, the estimated 
(M) are highly uncertain in many cases. Of particular im- 
portance is an estimate of the thermal response time-scale 
Tth of a neutron star after a significant change of (M). This 
can be easily obtained from Eq. JTJ after specifying Cv- 
Most of a neutron star specific heat is provided by its core 
and, for degenerate Fermions, CV oc T so we simply write 



2 This curve corresponds to the minimal cooling scenario, but 
any scenario fitting the isolated neutron stars data must have a 
similar L v . 



Cy = C-T erg KT 1 jYakovlev fc Pethicklliooi : IPage et al.1 
|2006| ). Then Eq. ^ gives 



d(T 2 
dt 



(7) 



where H and L are short term time averages of the heating 
term H and the cooling term L = L~, + L v . By "short term" 
time average we mean averaged over many accretion out- 
bursts but during a time span still much shorted that Tth- 
We assume that at time t = there is an abrupt change of 
H from Ho to Hi, and that at times t < the star was in a 
stationary state at a temperature To. Then at t < one had 
Lo = Hq. At t > the star begins to react and the resulting 
evolution of T 2 is illustrated by the thick curves in Fig. 
At times runs T will tend toward a new stationary state: 
T — ?> Ti and L —¥ L\ = Hi. One can estimate the thermal 
response time r t h by using a straight line (dotted line in the 
figure) from the initial position at t = and T = To till it 
reaches Ti at time Tth- The slope of this T 2 trajectory at 



t = is simply 2(Hi 
evolved. This gives us 



Lq)/C since L = Lo has not yet 



Tth 



C{T? 



T 2 : 



_ C{T 2 -T§) 
2(Hi-L ) 2{Hi-H ) 



were we also used that Lo = Ho- If Hi <C Ho or Hi ^> Ho 
this expression for Tth will be dominated by the high H , high 
T, term and we can write 



Tth 



2{if)high 



(8) 



Models show that C « 10 when T is measured in Kelvin, 
with only a small dependence on the mass of the neutron 
star. However it ca n be reduced by u p to a factor 10 in the 
presence of pairing (|Page et alj|2006h , i.e., C — > /sfC, with 
fsF ranging from 1 (no pairing) down to 0.1 (maximum 
extend of superfluidity/superconductivity). (H) is obtained 
from (M) with Eq. ([2]). Notice that since L 7 oc T 2 one has 
the fortuitous, but useful, result 



Tth 



2 x 10 4 fsF 



10 32 erg s" 



lQ- n M yr- 

(M) 



yrs. 



(9) 



Several tracks (the dotted lines) for various values of Tth are 
also shown in Fig. [T] 

Notice that the obtained Tth is just the Kelvin-Helmoltz 
time-scale t K h = Eth/L since E t h = J C v dT = ^CT 2 . 
However, when a star transit between a high and a low (M) 
states we have two tkh, and our analysis tells us that the 
relevant one is the high state one. If the energy loss is due to 
neutrinos the high state tkh is the shortest one since L„ oc 
T 8 or T 6 and thus tkh oc T" 6 or T" 4 . This is intuitively 
clear from Fig. [2] if the star starts from a low (M) its heating 
is driven by the new high (M) while if the star start from 
a high (M ) its cooling is driven by the remaining previous 
high L v . However, if the cooling is driven by the photon 
luminosity L 7 oc T 2 then Tth is essentially independent of 
<M): Eq © and © give r t 7 h « 10 5 f SF yrs, as is confirmed 
by the Tth line in Fig. [T] In both cases, it is simply the 
shortest tkh that determines the time-scale. 

Returning to the issue raised at the end of Section 12.11 
about the mass of these neutron stars, we see from Fig. [T] 
that the ones requiring fast neutrino emission, and hence a 
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Xth t 

Figure 2. Schematic evolution of star's core T 2 after an abrupt 
increase (top panel) or decrease (lower panel) in (M). See Sec- 
tion [2]2] for details. 

high mass, have a r t h ~ 10 3 — 10 4 yrs. Postulating that they 
have a low L q because (M) was much lower in the past and 
suddenly increased recently would require that this increase 
occurred in the last 10 3 — 10 4 for all of them. Such a claim 
is hardly sustainable and the most natural interpretation of 
the large spread in observed L q , for a given (M), is the result 
of a large range of neutron star masses: more massive stars 
should have a larger neutrino emission efficiency and hence 
appear fainter. 

2.3 Evolution of neutron stars with varying M 

When (M) changes, the L q of the neutron star will slowly 
evolve on a time-scale Tth- It will follow a trajectory, in the 
Lq vs (M) plane, given by Eq ((SJ) which is just one of the 
many curves plotted in Fig. [1] depending on its neutrino 
emission efficiency (i.e., and therefore very likely its mass). 

If neutrino cooling dominates the energy balance, i.e., 
L v = (H), we have for the star's core temperature T: T 2 oc 
(H) 1/3 for fast (T 6 ) neutrinos and T 2 oc (H) 1/4 for slow 
(T 8 ) neutrinos. Since L q oc T 2 we obtain for fast neutrino 
emission 

L q oc {M) 1/3 (10) 
and 

L q oc (M) 1/4 (11) 

for slow neutrino emission. For low enough (M) these neu- 
trino trajectories curve down and merge into the pho- 
ton cooling trajectory L q — (Qnuc/m u ){M} . These simple 
power-law behaviors are clearly seen in Fig. [1] However, for 
a very rapid change in M the neutron star may be found 
"off-trajectory" for a time ~ rth that can be as large as 10 5 
years for very faint systems. 

Notice that from the theory point of view the photon 
trajectory is the cleanest one: it only depends on Q nuc , 
Eq. Q, which is known within a factor two, 1-2 MeV, 
within the deep crustal heating model. The neutrino trajec- 



tories (which depend on the core temperature T) are plagued 
with the uncertainty on how to relate T e with the core T. 
This T e — T relationship depends on the outer layers chem- 
ical composition and the star's surface gravity: for a given 
neutron star model (encapsulated in our parameters N and 
C for the neutrino luminosity and the specific heat) the 
predicted L q can vary by almost on e order of magnitude 
because of the T e — T relationship (|Potekhin et al.l 1 19971 : 
lYakovlev et~"ai1l2004l ). 

2.4 Testing the deep crustal heating model 

So far this deep-crustal heating model has mostly been 
tested using sources within a relatively small (M) range, 
between 5 x 1CT 12 to 5 x 1CT 10 M yr _1 , as seen from 
Fig. [T] For higher (M) the model will be difficult to test 
because it straddles in the range of the persistently bright 
sources and quiescent measurements are impossible for such 
sources. However, at the low (M), between 1 x 10~ 13 erg s _1 
and 5 x 10 12 erg s _1 , it is likely that sources can be added; 
in particular the neutron-star VFXTs can have (M) in this 
range (see section [3TTJ . We now discuss the prospect of us- 
ing those systems to test th e deep-crustal h eating model (see 
also the brief discussion in IWiinanddl2008h . 



3 THE QUIESCENT LUMINOSITY OF 
NEUTRON-STAR VFXTS 

From the results of the previous section, illustrated in Fig.[T] 
we can estimate the expected quiescent luminosities L q of 
VFXTs from estimates of their long term average accre- 
tion rates. However, it is notoriously difficult to estimate 
(M) for any X-ray transient and even more so for VFXTs 
because they are diffi cult (despite the sens itivity of X-ray 
instruments in orbit; IWiinands et al.l [2006) to detect and 
therefore most of their outbursts are likely missed. For the 
deep-crustal heating model, the (M) which is needed is the 
M averaged over a time Tth that (as can be seen from Fig. [1} 
ranges between 10 4 and 10 5 years at M ~ 10~ 12 Mq yr" 1 
and L q ~ 10 31 erg s" 1 . 

3.1 Estimation of (M) in VXFTs 

An estimation of an upper limit on (M) for the VFXTs 
can be obtained using some simple assumptions. As already 
stated in the introduction, VFXTs have peak X-ray lumi- 
nosities (2-10 keV) < 1 x 10 36 erg s" 1 . Thi s limit was 
chosen in the classification of IWiinands et al.l (|2006h be- 
cause it roughly corresponds to the sensitivity limit of the 
past and present X-ray all-sky monitors in orbit (e.g., Bep- 
poSAX/WFC, RXTE/ASM, Swift/B AT, Integral, MAXT) 
for sources at 8 kpc. Typically those instruments are sen- 
sitive to X-ray outbursts which have a peak flux (2-10 keV) 
above 10 mCratlf], which corresponds to a 2-10 keV flux of 
~ 2 x 10 -10 erg s _1 cm -2 and a limiting 2-10 keV luminosity 

3 Some instruments or programs are more sensitive (e.g., 
the RXTE/PCA bulge sc an project reached about 1 mCrab; 
ISwank fc Markwardtll200lt) . However, for the purpose of this pa- 
per, we assume a conservative limit of 10 mCrab. 



© 0000 RAS, MNRAS OOP, HUTU 



6 Wijnands, Degenaar, & Page 



sensitivity of 1 x 10 36 erg s _1 (for 8 kpc; which was assumed 
the typical distance towards VFXTs since most have been 
found in the Galactic bulge). 

Typica lly the bolometric lum inosity is then a factor 2-3 
higher (see lin 't Zand et al.ll2007T l and thus the upper limit 
on the bolometric peak luminosity would be < 3 x 10 36 erg 
s _1 . Assuming perfect efficiency of the accretion process, 
this accretion luminosity is related to the accretion rate in 
outburst M through L — GJ1 ^ M , with G the gravitational 
constant, M the mass of the neutron star and R the radius 
of the star. Here we use a "canonical" neutron star with a 
mass of 1.4 solar masses and a radius of 10 km. The lu- 
minosity upper limit then result in an upper limit on the 
peak mass accretion rate during outburst of < 3 x 10~ 10 
Mq yr -1 . Assuming that the duty cycle DC (with DC = 
to I {to +t q ), t being the outburst duration time and t q the 
quiescent duration time) of VFXTs is between 1%-10% (as is 
typically observed for the recurrent bright transients), then 
(M) < 1(T 12 - 1CT 11 M© yr" 1 . This is a rather conservative 
upper limit because it assumes a "step-function" outburst in 
which the sources are always accreting just below the limits 
set by the all-sky instruments when they are in outburst. 
However many sys tems will have peak luminosi ties well be- 
low this level fe.g.. |Pegenaar fc Wi jnands 20l"(ih and during 
the outbursts the peak luminosities are only reached during 
a small fraction of the outburst, sim ilar to the outbur st pro- 
files of brighter transients (e.g., see lChen et al.|[l997l ). This 
will significantly lower this upper limit for the VFXTs. How- 
ever, this limit already shows that indeed VFXTs have very 
low (M); i.e., in the range which has hardly been used to 
test the deep-crustal heating model. 

Moreover, for several VF XTs more stringent cons t raint s 
can be obtained on their (MV lDegenaar fc Wijnands! (|2010l ) 
have estimated the (M) of several VFXTs near the Galactic 
center using 4 years of Swift/XKI monitoring data on Sgr 
A*. The majority of the sources considered in that paper 
had (M) < 2 x KT 12 M yr" 1 and were typically in the 
range of 1CT 13 - 10" 12 M© yr" 1 . This demonstrates that 
at least some VFXTs have extremely low (M) and it might 
be possible that a large fraction of the VFXTs have similar 
(M). 



3.2 Estimation of L q in VFXTs 

A naive estimation of the possible range of L q 's for neutron- 
star VFXTs would simply consider the various trajectories 
of Fig. [T] we can consider the sample of objects depicted in 
Fig. [I] and move them along their respective cooling trajec- 
tories into the range of the (very low) (M) above estimated 
for the VFXTs. There are, however, two immediate issues 
with this extrapolation: 

1) How representative are the above estimated (M) of 
VFXTs with their real long term (M)? 

2) Since the spread in the L q observed for the brighter 
transients exhibited in Fig. [T] is likely due to a large spread 
in neutron star masses (as argued above) how representative 
is the mass distribution of the bright transients for the mass 
range of the VFXTs? 



3.2.1 The (Af) versus Tth issue 

The above limits on the (Af) assumes that the current ob- 
served behavior of VFXTs is representative of their gen- 
eral behavior. However, this does not need to be true, and 
it might be that some fraction of the VFXTs are actu- 
ally sources which usually accrete at much higher rates 
and we only observe them during a small period in which 
their (M) is much less (although see the d iscussions in 
iKing fc Wijnands! 120061 ; IWiinands et"aH l200rj . demonstrat- 
ing that this is likely not true for the majority of the 
VFXTs). 

The issue is how long has this period of very low 
(Af) lasted compared to the neutron star relaxation time Tth- 
As discussed in Section [2~3l if the time-scale of the evolution 
of (Af ) is larger than the corresponding Tth the neutron star 
will evolve along one of the trajectories exhibited in Fig. [1] 
this is the Case B discussed below. However, in the opposite 
case we have: 

Case A: The VFXTs in this class are normally bright 
transients or even bright persistent sources and we only 
observe them during a short-lived VFXT episode, shorter 
than their Tth- This is not unfeasible, because several bright 



well (e.g.. iLinares et al.l 


20071: iDegenaar & Wijnands! 120091: 


iFridriksson et al.l 2011: 


Degenaar et al.l 12012a). In other 



words, the estimated (Af ) is rather an Af in the notations 
of Section \2. 21 and the star has not yet adjusted to the new 
low Af : the star's present L q is still determined by its previ- 
ous higher (Af ) and is "off-trajectory" . In the case the neu- 
tron star is not too massive and was previously on a slow 
neutrino cooling track its present L q can be in the range 
of 10 33-34 erg s _1 : such high L q can be sustained for more 
than 10 yrs, the corresponding Tth of bright transients on 
the slow cooling trajectories {rth being determined by the 
value in the high (Af) regime). With such an L q and an 
(Af ) in the VFXT regime it would be located in the diagram 
of Fig. [T] above the photon cooling line. This is a region that 
is unaccessible for a star in thermal equilibrium. 

Case B: The observed behavior of the VFXTs in this 
class is representative for their behavior on a long enough 
time that they are in a steady state, i.e., they can be cor- 
rectly located on one of the trajectories of Fig. [T] However, 
for the Case B systems, one has to consider the second issue, 
which we will address in the next section. 



3.2.2 The mass distribution issue 

Extrapolating the observed distribution of L q of bright tran- 
sients to the low (Af) regime could allows us to predict the 
expected distribution of L q for the VFXTs. However, iso- 
lated young neutron stars, that likely have a mass distribu- 
tion between 1.2 to 1.6 Mq, with an observed thermal lumi- 
nosity are well describe d by slow neutrino emission processes 
|Page et al.ll2004l2009l ). This is in sharp contrast with what 
is observed for the bright neutron-star transients. Most of 
them have an L q that requires fast neutrino processes, im- 
plying a mass distribution strongly skewed toward higher 
masses due to the long term accretion. As a consequence, 
the L q distribution of the VFXTs that have been in the 
VFXT phase for long enough that they are in a steady state 
(Case B) should show a significant imprint of their very long 
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term past accretion history and Case B can be divided into 
two distinct sub-cases: 

Case Bl: Some systems might have been a VFXT for 
a period longer than rth but before that they were similar 
to the brighter transients. Although originally the tempera- 
ture of the core should be relatively high, the time spend in 
the VFXT phase is long enough that the neutron star has 
adjusted its core temperature (i.e., it has become colder) to 
the very low accretion rate, but during the time spend before 
this phase (when the accretion rate was much higher) the 
source could have accreted enough material for it to cross 
the mass threshold beyond which fast core cooling could oc- 
cur. Therefore, those systems should be on one of the "fast 
neutrino cooling" trajectories of Fig. [T] and be significantly 
fainter in quiescence than predicted using the standard deep- 
crustal heating model, Eq. (Q. 

Case B2: Some systems might have been VFXTs 
throughout their life and the estimated (M) is representative 
for the (M) thr oughout their life (we ca ll them primordial 
VFXTs; see also lKing fe Wiinandj|200fj| . who talked about 
the possible existence of such systems). Therefore, since it 
is expected that LMXBs live for 10 8-9 years, only about 
10" 5 to 10~ 2 solar masses (assuming (M) = 1CT 13 - KT 11 
Mq yr -1 ; Section l3.ip has been accreted by the neutron 
stars. Therefore, it is quite likely that most, if not all, of 
these sources have not been accreting enough matter dur- 
ing their life to increase their masses (assuming they are all 
born with a mass of ~1.4 M©) above the threshold so that 
enhanced core cooling would become active. If true these 
systems should be on the slow neutrino cooling tracks of 
Fig. [T]that, in the estimated (M) range of VFXTs, merge 
with the photon cooling track. Therefore, they should agglu- 
tinate on the narrow "photon cooling" region and indeed be 
detectable at io 31-33 erg s _1 . However, this assumes that 
the standard deep-crustal heating is active in those sources, 
but in the next section we argue that this might not be the 



3.3 Non-standard heating in neutron-star VFXTs 

In the deep-crustal heating scenario, it is explicitly assumed 
that the origi nal, catalyzed crust is fu l ly rep l aced by an 
accreted one jHaensel fe Zdunikl [l990al . 1200a . l2008r i. The 
low-Z elements present (at high densities) in an accreted 
crust facilitate a significant amount of energy release due to 
the pycnonuclear reactions. Those pycnonuclear reactions 
occur in the density range from 10 12 up t o ~ 5 x 10 13 g 
cm" 3 jHaensel fc Zdunikl 1990 j , 12001 moj ). It is quite pos- 
sible that among the primordial VFXTs (the above Case B2 
sources) a group of systems exist in which the crust is not 
replaced to this depth and a partly accreted and a partly 
original (albeit further compressed) crust (a hybrid crust) 
might be present. 

To make this last assertion more quantitative, we show 
in Fig. [3] estimates of the amount of matter that needs to 
have be accreted in order to replace the original catalyzed 
matter. The needed accreted mass is simply 
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where y is the accreted column density. This column den- 
sity determines the pressure reached by the accreted mat- 
ter since, from hydrostatic equilibrium, P — gy where 



Figure 3. Estimate of the density reached by accreted mat- 
ter as a function of total accreted mass AMacc for 4 different 
neutron star masses and 3 different radii. Lines show the values 
obtained assuming a catalyzed crust EOS while the background 
shaded areas show the corresponding range of values for a wholly 
accreted crust. The vertical "pycnonuclear" double-arrow marks 
the density range in which pycnonuclear fusions are expected to 
take place, and significantly heat the neutron star, once catalyzed 
matter has been replaced by accreted matter. See Section 13.31 for 
details. 



g = e GM /R is the gravity acceleration in the crust and 
e* = (1 - 2GM/Rc 2 )~ 1/2 the red-shift (G being the gravi- 
tational constant and c the speed of light). For a given pres- 
sure we can obtain the corresponding density p using a crust 
equation of state (EOS). The EOS of a hybrid crust remains 
to be calculated but we can bracket it between the EOS of a 
catalyzed crust and a wh olly accreted o ne. F or a catalyzed 
crust we use the EOSs of | Haensel et 3 (| 19891 1 for the outer 
crust and lNegele fc Vautherinl ( 1973T ) for the inner crust. For 
an ac creted crust we use the model of lHaensel fc Zdunikl 
(|2008r j. The results of Fig. [3] show that, fortunately, the 
needed AM ncc does not depends strongly on the assumed 
EOS. Pycnonuclear reactions provide about 60-70% of the 
total Q„uc and one sees that, for a low mass (1.2 Mq) ex- 
tended (14 km radius) neutron star, about 2 x 10 -2 Mq of 
matter needs to be accreted so that Q nU c reach its optimal 
value while about 2 x 10 -4 Mq is needed for the accreted 
matter to reach the threshold of the first possible reaction. 
For a heavy (1.8 Mq) compact (10 km radius) neutron star 
the corresponding numbers are approximately 1 order of 
magnitude smaller but such a heavy star will likely undergo 
fast neutrino cooling and have such a low L q in the (M) 
regime of primordial VFXTs that it will be practically un- 
observable. We note that if such heavy neutron stars could 
indeed be present in some primordial VFXTs that those 
neutron stars have to have been born this massive. 

The effect on the heating is not clear. It depends on 
to what depth the crust has been replaced since pycnonu- 
clear reactions occur down to a density of ~ 5 x 10 13 g cm -3 



but can start already at significantly lower density (see, e.g., 
iHorowitz et al.|[2008f ). In particular most of the heat due to 



© 0000 RAS, MNRAS OOO.HTfMl 



8 Wijnands, Degenaar, & Page 



pycnonuclear reactions is released in the density range of 
10 12 to 10 13 g cm -3 (e.g. jHaensel fe ZduniklbOOct ) and sig- 
nificantly less mass has to be accreted to replace the crust 
to those densities (Fig. |3). However, in the most extreme 
case in which the total amount of accreted matter is only 
10~ 5 Mq then only the outer crust, with densities below 
the neutron drip, has been fully replaced, significantly in- 
hibiting pycnonuclear reactions to occur in the inner crust. 
Therefore, there might be a sub-group among the primordial 
VFXTs for which the heating is significantly reduced which 
would make them even fainter than already inferred from 
their low (M) . 



4 DISCUSSION 

We have estimated the quiescent thermal luminosity of 
neutron-star VFXTs in order to determine if they can be 
used to test the deep-crustal heating model in an hardly ex- 
plored (M) regime. Unfortunately a conclusive answer can- 
not be give due to the large uncertainties in our knowledge of 
the accretion rate history of VFXTs. The (M) of the source 
during the last several thousands to tens of thousand years 
determines how much heat has be deposited in the neutron 
star over that period and therefore the thermal state of the 
star. However the long-term history over the lifetime of the 
binary determines the amount of matter accreted and there- 
fore if enough matter has been accreted to trigger enhanced 
neutrino emission processes in the core and if enough mat- 
ter is accreted to allow the activation of all pycnonuclear 
heating reactions in the inner crust. 

This last point arises because it is well possible that 
the amount of matter which primordial VFXTs have ac- 
creted during their lifetime is not enough to fully replace 
the original crust, leaving a crust which is partly replaced 
by accreted matter and partly still contains the original, al- 
beit further compressed material. It is unclear how such a 
hybrid crust would react to the accretion of matter and how 
this would effect the thermal state of the neutron star. Likely 
less heat is produced because not all pycnonuclear reactions 
can occur, but it is not clear if other properties of a hy- 
brid crust are also significantly different compared to a fully 
accreted crusts, such as the thermal conductivity. Besides 
obtaining more observational data to constrain the models, 
detailed theoretical calculations have to be performed to in- 
vestigate the heating and cooling in neutron stars which have 
a hybrid crust. In particular it is important to investigate 
different crustal compositions with a variety of amount of 
m atter accret ed (e.g., an update study of the one performed 
by |Satdll979T ). This problem might not only be interesting 
for VFXTs, but also to other types of neutron stars because 
VFXTs might not be the only sources which harbor neu- 
tron stars with hybrid crusts (see Section ^. ip . Furthermore, 
VFXTs might be important to understand low-luminosity 
X-ray sources in globular clusters (Section 14. 2[) . In addition, 
they might form an interesting group of sources to try to 
study cooling of the neutron star crust (Section l4.3|) after it 
has been heated during outbursts. 



4.1 Additional potential sources without fully 
accreted crusts 

Despite that it is generally accepted that most neutron-star 
LMXBs are rather old systems with ages of 10 8-9 year, 
there are individual sources which likely are much younger. 
One relatively young system might be the recently discov- 
ered transiently accreting 11 Hz X-ray pulsar IGR J17480- 
2446 i n the globular cluster Terzan 5 (also called Terza n 
5 X-2; IStrohmaver fc Markwardt|[201ol ; iPapitto et alj|201lf ). 
This system is an unusual LMXB because it was expected 
that the neutron stars in LMXBs should have spin periods 
< 10 millisecond because t hey are spun up by the accretion 
of m atter (see review by iBhattacharva fc van den HeuveJ 
Il99ll) . The slow spin period of IGR J17480-2446 is enig- 
matic and it has been hypothesized that this is due to the 
fact that the system has so far only spend a relatively brief 
time in the Roche- lobe overflow phase (10 7 to 10 8 years; 
iPatruno et al.|[2012r >. If true, this system might be an exam- 
ple of systems which do not have a fully replaced neutron 
star crust. 

The mass accretion rate of this source during out- 
burst has been estimated to be 3 x 1CP 9 Mq yr _1 
l|Degenaar fe WiinanddfeoilaT l. The duty cycle of this sys- 
tem is poorly constrained but if we assume again values 
of 1%-10% we obtain a time-averaged accretion rate of 
3 x 1CT 11 - 3 x 10~ 10 M yr" 1 . Combined with the expected 
age of the accretion phase this results in a mass accreted on 
the neutron star of 3 x 1CT 4 to 3 x 1(T 2 M Q . Although 
the maximum amount of matter accreted would indicate 
that the full crust is replaced, it is also quite possible that 
the neutron star in th i s syste m has a hybrid crust as well. 
iDegenaar fc Wijnands! (1201181 ) found the quiescent counter- 
part for this source to be rather cold, significantly colder 
than expected using standard heating and cooling theory. 
They suggested that in its neutron star enhanced core cool- 
ing processes might be active although, as also shown above, 
probably not enough matter has accreted on the star for the 
star to have become massive enough to allow such processes 
to occur in the core. Alternatively, they suggested that the 
duty cycle might be extremely low, of the order 0.1%. 

Although not impossible, this duty cycle seems very 
low (and possibly improbable in the disk instability model) 
and therefore we suggest an another possible reason why 
the source is so faint in quiescenc^S due to the presence of 
a hybrid crust, not all the heating reactions can occur in 
the crust and therefore less heat has been deposited in the 
neutron star to heat it up to the expected temperature as 
inferred from its (M). This conclusion still holds when also 
taking into account that before th e Roche-lob overflow phase 
a wind-accretion phase occurred. IPatruno et al.l (|2012T ) esti- 
mated that the mass accretion rate in that phase would at 
most be 10~ 13 - 10~ 12 M yr~\ This phase could have 



4 We note that the thermal quiescent luminosity of the source is 
still well within the range observed from other quiescent neutron- 
star LMXBs which would suggest that the source is not special. 
This could indicate that the same physical processes are at work 
in this source as well as in the other sources. This would also 
satisfy the principle of Occam's razor, by not having to have to 
postulate several mechanisms why certain quiescent LMXBs arc 
colder than expected by the standard theory. 
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lasted 10 7 " 8 year and thus at most 10 6 - 1CT 4 M could 
have been accreted. 

If this is the correct explanation for why the neu- 
tron star in IGR J17480-2446 is colder than expected, one 
has to wonder if a similar argument might also hold for 
other systems which have found to be too cold. For ex- 
ample, the neut ron star in SAX J1808 . 4-3658 seems to be 
extrem ely cold jCampana et al.ll2002al ; iHeinke et all 12001 
l2009al). Its (M) has bee n estimated to be ~ 10" 11 M 
yr _1 (|Heinke et alj [2007) and if it has lived shorter than 
10 s years the neutron star should have a hybrid crust. How- 
ever, this system is an a ccreting millisecond pulsar wit h 
a spin period of 401 Hz jWiinands fc van der Klisl ll99ST ) . 
This means that a significant amount of matter has to have 
been accreted by the neutron star to spin it up to this 
spin freque ncy. Typically, the cal culations show that at up 
to 0.1 Mq (|van den Heuvell 1987) is needed to accomplis h 
this (see review by iBhattacharva k, van den Heuvellll99j ). 
Therefore, in SAX J1808. 4-3658 the neutron-star crust will 
have been fully replaced, which strongly indicates that in 
the past the accretion rate of this system was considerably 
larger than its current inferred (M). Another source which 
might be relatively young is Circinus X-l. The age of this 
system is not known, but it has been suggested to be rather 
young (of the order of < 10 4 ~ 5 years; see the discussion in 
IClarkson et al.l [20041 ). Despite that it can accrete on occa- 
sions at very high accretions rates (up to > 10 M© yr _1 ), 
this age (if confirmed) is sufficiently low that very likely not 
the complete crust has been replaced. If the source would go 
fully quiescence, it would be very interesting to determine 
the quiescent luminosity of the neutron star in this system. 



4-1.1 Neutron stars in high-mass X-ray binaries 

In high-mass X-ray binaries (HMXBs) the neutron star ac- 
cretes either from the strong stellar wind of the companion 
(e.g., a supergiant star) or from the decretion disk of a B 
type star, which is typically observe d to be of t ype B0-B2 in 
Be/X-ray transients (see review by |Reig||201ll ). Such early- 
type B stars only live between 10 to 30 million years. Typi- 
cally in Be/X-ray transients the sources have outbursts with 
X-ray luminosities of io 36-37 erg s _1 (corresponding to an 
outburst accretion rate of 10~ 10 to 10 -9 Mq yr _1 ) when the 
neutron star moves through the decretion disk at periastron 
passage. For some sources this occurs once every orbital pe- 
riod (resulting in periodic outbursts; called type-I outburst) 
but other sources are only occasionally in outburst. There- 
fore, it is unclear what fraction of the time the neutron star 
is actually accreting, but when assuming again a duty cycle 
of 1%-10%, this would result in a (M) of 10~ 12 to 10~ 10 M Q 
yr -1 and a total amount of mass accretion throughout the 
life time of the system (assuming the system was a Be/X-ray 
transient for the whole life of the B star which might be a 
significantly overestimation of the duration of this phase) of 
10~ 5 to 3 x 10~ 3 Mq. Thus, it is quite possible that also the 
neutron stars in some Be/X-ray transients have a hybrid 
crust. We note that some systems also exhibited so-called 
type-II outbursts which are much brighter (peak luminosi- 
ties of 10 38 erg s _1 ) which can last for weeks to months but 
they are very infrequent and not all systems exhibit them. 
Therefore, we do not expect that those type of outbursts 
will affect our main conclusion significantly. 



Another class of HMXBs transients are the supergi ant 
fast X-ray transients (or SFXTs; see, e.g.. ISidolil l201ll ) in 
which the neutron star transiently accretes from the vari- 
able dense wind of a supergiant star. However, also very 
likely in those systems the neutron star has only a partly 
replaced crust because the supergiants only live very short 
and despite that the outbursts of those systems can be very 
bright (10 38 erg s _1 ) they are very brief, very infrequent, 
and most of the time the neutron star is only accreting at 
much lower rates or not at all. Although, if before the su- 
pergiant phase the neutron star was also already accreting 
significan tly from the com panion star (e.g., during an earlier 
Be phase; iLiu et alllioTll ). then more of the original crust is 
replaced. 

For the neutron stars in HMXB transients (M) is typ- 
ically higher than inferred for the VFXTs (they are typi- 
cally more in the range observed for the ordinary LMXB 
transients). Therefore, it is expected that if standard heat- 
ing and cooling o ccurs in those systems (as suggested by 
I Brown et al. I [19981 ). that their thermal emission should be 
readily detectable in quiescent. Enhanced core cooling is not 
expected because they should be relatively light weight neu- 
tron stars since little matter has been accreted (although it 
might be possible that some systems are b orn with massive 
neutron stars; see, e. g., iBarziv et al1l200ll) . In contrast, the 
heating might be affected by what type of crust is present 
(i.e., fully accreted or hybrid crust) and HXMB transients 
might be very good candidates to investigate the effect of 
hybrid crust on the thermal properties of the neutron star. 
However, the situation for those sources might be compli- 
cated by the much stronger magnetic field in those systems 
(10 12-13 Gauss) compared to those of the neutron stars in 
LMXBs (10 8-9 Gauss). It is unclear how strong the effects 
of these stronger magnetic fields are on the heating and 
cooling of the neutron stars and other related properties 
(e.g., the thermal conductivity which is severally affects by 
super strong magnetic fields of > 1 13 Gauss and there- 
fore likely also by sl ightly lower fields; iPotekhin et al. I ll999l ; 
lAguilera et al.ll2008l ). More detailed theoretical calculations 
have to be performed to determine the effect of the mag- 
netic field, in combination with the exact composition and 
structure of the (possible hybrid) crust. 

Observing HMXB transients in their quiescent state 
could be very useful in this aspect. However, the number 
of neutron-star Be/X-ray transients so far s tudied in quies- 
cence is rather limited (fo r a source list see iRutledge et all 
120071 ; iTomsick et all [20111 ). So far, the obtained picture is 
complex. Some systems (like, e.g., EXO 2030+375) always 
remain rather bright in-between outbursts (> 10 35 erg s , 
basically they never transit to quiescence). However, the 
majority of systems h ave quiescent luminosities between 
10 32 and 10 34 erg s" 1 (|Rutledge et al.ll2007l ; ITomsick et all 
l201ll ). Spectral analysis demonstrates that some systems are 
still very hard in quiescence with power-law indices near 1 
or even lower (simil ar to what often is seen in outburst; 
IRutledge et al.|[2007i . e.g., ), while others are soft er with in- 
dices even up to 2.6 (e.g.. lcampana et al.|[2002bl ). 

Although the quiescent data are usually not of very 
high quality, several sources do not show pulsations in quies- 
cence which might indicate that indeed the accretion down 
to the surface has halted in those systems (|Campana et al.l 
l2002bl ; IWilson et ai]|2005l ) . However, in a few other systems 
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pulsations could still be detected in quiescence demonstrat- 
ing that in those systems either some of the matter still 
reaches the neutron-star surfaces or the pulsations are in 
some way caused by the interaction between the magnetic 
field (which is rotating with the neutron star ) and the accre- 
tion of matter down to the m agnetosphere (Rutlcdg c et al.l 
l2007l : lMukheriee fe Paulll2005l ). 

For some syste ms it has been suggested 
jCampana et all l2002bl l that indeed the emission we 
observe is due to the cooling of the neutron star and not 
due to some sort of accretion process, however, the evidence 
is not conclusive due to the statistical quality of the 
data. Furthermore, the possibility that they might harbor 
a neutron sta r with a hybrid crust was not considered. 
iTomsick et"aH (|201lh discussed possible reason why the 
candidate Be/X-ray transient IGR J01363+6610 could not 
be detected with Chandra in its quiescent state (e.g., the 
system containing a black hole instead of a neutron star), 
but in light of the above discussion we suggest that the 
possibility should be considered that this source might still 
harbor a neutron star but one with a hybrid crust which 
inhibits significant heating of the neutron star. 

The situation for SFXTs is similar to that of the Be/X- 
ray transients with only a handful of SFXTs studied in qui- 
escence. Also tho se systems show a variety in quiescent be - 
havior (see, e.g.. lin't Zand! 120051 ; IBozzo et alJboid . I2012T ). 
A systematic and homogenous study of many more HMXB 
transients (both SFXTs and Be/X-ray transients) in quies- 
cence is needed to understand fully how they can be used to 
study the deep-crustal heating model. A survey of 16 con- 
firmed neutron-star Be/X-ray transients in their quiescent 
state has recently been accepted (PI: Wijnands) which will 
give more insight into this issue. 



4.2 VFXTs in globular clusters 

Many faint X-ray sources have been found in the Galactic 
globular clusters, and a large number are li kely associated 
with neutron-star X-ray transients (see, e.g.. IVerbunt et all 
1 19841 ). But the lack of a significant number of outbursts 
from those sources has lead to suggestions that maybe those 
sources are associated with VFX Ts whose outburs ts where 
missed by the all-sky monitors (|Wiinandsl |2008| ). As es- 



timated in section 13.21 the quiescent X-ray luminosity of 
VFXTs in the standard deep-crustal heating model would 
be in the range io 31-33 erg s _1 and indeed, if the stan- 
dard heating and cooling processes occur, a large fraction 
of the candidate quiescent LMXBs could be associated with 
VFXTs. However, as also explained in section 13.21 and 13.31 
the quiescent luminosity of VFXTs could be significantly 
lower than expected in the standard model and therefore it 
is unclear if this conclusion still holds. Moreover, even in the 
standard model it is unlikely that the VFXTs are associated 
with the candidate quiescent LMXBs in globular clusters. 

To demonstrate this, we rewrite the time-averaged ac- 
cretion rate into 



(M) 



(M )t + (M q )t q _ ^ ^ 



to 



(M a ) DC (13) 



to ~t~ tq to ~\~ tq 

with (M ) the time-averaged accretion rate in outburst and 



{M q ) the time-averaged accretion rate in quiescence. Equa- 
tion ll3l assumes that {M )t ^> {M q )t q , which is usually true 
but might not if t q S> t , thus for systems with a extremely 
low duty cycle. Using equations [2] and 1131 and assuming 



(H) (thus standard, slow cooling) one obtains 



DC = 



10" 



M 



O yr 



1.5 MeV 



10 33 erg s _1 



> 



{Mo) 



Qnuc 

1.5 MeV 



3 x 10 34 erg s" 1 



(14) 
(15) 



which assumes (M a ) < 3 x 10~ 10 Mq yr" 1 (the limit set by 
the all-sky monitors). 

Typically quiescent LMXB candidates in globular clus- 
ters have bolometric X-ray luminosities of 10 32 ~ 33 erg s - 1 



with the brightest being ~ 3.6 x 10 33 erg s 1 (|Heinke et al.l 
120031 , about a third of the sources listed in that paper 



have an X-ray luminosity > 1 x 10 erg s ). This results 
(using equation I15|l in limits on the duty cycle of DC > 



0.0025 - 0.005 for L a = 1 x 10 32 erg s - \ DC > 0.025 



0.05 for L a 



1 x 10 33 erg s" 1 and DC > 0.1 



0.2 for 



that we have assumed that Qnuc must be between 1 and 
2. Putting reliable observational constraints on the duty cy- 
cle of possible VFXTs in globular clusters is difficult. How- 
ever, currently there are about 30 quiescent LMXB candi- 
dat es identified for whic h no outburst^] have been seen yet 
(see lHeinke et all (|2003l) for a list, seve ral additional source s 
have found si n ce that publ ic ation; e.g.. [ Lugger et al.l (|200?f ) ; 
iGuillot et al] (|2009l . l201ll ); iMaxwell et al.l i|2012l )). During 
the Chandra and XMM-Newton observations which detected 
those sources, they were not in outburst. Therefore, assum- 
ing that all sources are from the same population of tran- 
sient typ e, the duty cy cle of those systems is < 1/30 = 0.033 
(see also lHeinkell2010l ). We note that a number of clusters 
have multiple Chandra and XMM-Newton observations dur- 
ing which those sources were not seen in outbursts. However, 
those additional observations do not constrain the duty cy- 
cle further because it is quite possible that the quiescent 
duration is significantly longer than the sampling time scale 
meaning that the observations probe the same quiescent pe- 
riod and therefore are not statistically independent (basi- 
cally many outburst-quiescent cycles must have passed for 
the sampling, if performed with random time-intervals be- 
tween the observations, to become independent). We note 
that the uncertainties on the duty cycle inferred from ob- 
servations is likely to be large, but the duty cycle has to 
be >10% to explain the brightest sources (in the standard 
deep-crustal heating model) which seems unlikely. Further- 
more, similar low dut y cycles were inferred for VF XTs near 
the Galactic center ( Degenaar fc Wijna nds 2010), so it is 
quite possible that they indeed have such low duty cycles. 
The conclusion from the above exercise is that VFXTs 



5 This excludes the two sources (IGR J17480-2446, Swift 
J174805. 3-244637) in Terzan 5 which were previo usly identified as 
candidate quiescent LMXBs llHeinke et al.|[2~006l) but which have 
now been shown to be associated with bright transients. How- 
ever, the final conclusions does not depend very sensitive on how 
many quiescent LMXBs are currently known but have not been 
associated yet with outbursts. 
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can still be associated with some of the quiescent LMXB 
candidates, but only with the faintest sub-set of the group 
and likely not with the brightest objects (> 10 33 erg s" 1 ; 
especially not those of a few times 10 33 erg s _1 ; note that 
this assumed standard heating and cooling). This is sup- 
ported by the fact that the three VFXTs currently known 
in gl obular clusters (M15 X-3 NGC 6440 X-2, a nd NGC 
6388; lHeinke et alj|2009bl , bold ; [Bozzo et alj|201lh all have 
quiescent luminosities ( well) below 10 32 erg s _1 (only M15 
X-3 has been detected; iHeinke et al.|[2009bl ). It seems that 
such faint transients are indeed present in globular clus- 
ters, but that their quiescent luminosities is really low. It is 
worth noting that those are the only systems for which use- 
ful constraints have been set on the quiescent properties of 
neutron-star VFXTs. The Galactic disk sources are usually 
too much absorbed (see also Section ^. 3[) for any useful con- 
straints indicating that globular clusters are prime targets to 
study quiescent VFXTs. Sensitive monitoring programs (see 
IWiinands et al.ll2012i; lAltamirano et al.ll2012f ) are needed to 
detect the very- faint outbursts of VFXTs in globular clusters 
and combined with rapid fol low-up observations using Chan- 
dra (|Altamirano et al1l2012r i to determine the exact position, 
the quiescent counterpart can then be studied in archival 
Chandra data or in newly proposed Chandra observations. 

It might be possible that we have underestimated the 
amount of energy released per accreted nucleon and there- 
fore the systems should be more luminous per accreted nu- 
cleon than assumed (e.g., evidence has been reported that 
in some systems extra heat in shallower layers in the crust 
must be li berated to explain the qui e scent properties of thos e 
systems; iBrown fc Cummindl2009l ; iDegenaar et al1l2011rj ). 
However, this should also be true for the other quiescent 
LMXBs, which do not require more energy per accreted nu- 
cleon to explain the base quiescence level and most are ac- 
tually too cold to be explained using the standard model. 
In addition, there is no theoretical basis for a significantly 
larger amount of energy released per accreted nucleon and 
therefore we consider this possibility unlikely. 

Therefore, we conclude that there are indeed VFXTs 
in globular clusters, but that they can likely not be as- 
sociated with the brightest amongst the quiescent LMXB 
candidates (> 10 33 erg s" 1 ) and likely also not with the 
slightly fainter ones (10 32 ~ 33 erg s" 1 ). It is more likely 
that those quiescent LMXB candidates are associate with 
bright transients similar to the transients IGR J17480- 
2446 and Swift J174805. 3-244637 in Terzan 5 which were 
as well identified pr eviously as quiescent LMXB candidates 
llHeinke et al.ll2006T) before they exh i bited their bright out - 



bursts jDegenaar fc Wiinandsll2011al ; IWiinands et al.ll2012t) . 

However, also those transients should have low duty cycles 
in order to be consistent with the observations of the lack of 
outbursts (either very-faint or bright). More and more evi- 
dence becomes available that there exist a class of transients 
which might have indeed very low duty cycles. This might 
not be unexpected because their is a strong selection effect in 
favor for discovering new transients with a high duty cycle. 

4.3 Final remarks 

The above discussion has demonstrated that VFXTs could 
be extremely faint in quiescence but some of them could 
reach quiescent X-ray luminosities of 10 32 to even 10 33 



erg s _1 depending on their accretion rate history and the 
amount of matter accreted. So, it would be interesting to ob- 
serve quiescent VFXTs and determine their quiescent prop- 
erties. Sadly, two main factors hamper the study of the ther- 
mal cooling emission of quiescent neutron-star VFXTs. First 
of all, many quiescent LMXBs exhibit besides the thermal 
component a hard, non-thermal component above 2 keV. 
The origin is not well understood but its presence inhibits 
the most accurate study of the thermal component and in 
some systems only this hard component can be detected. 
In particular, there seems a trend which indicates that the 
fainter a source is in quiescence, the larger the contribution 
is of this n on-thermal compon ent to the 0.5-10 keV X-ray 
luminosity |jonker et al.|[2004l ) . If this trend is also valid for 
the quiescent state of VFXT, then this will make it difficult 
to maybe even impossible to put significant constraints on 
the thermal component in those systems. 

In addition, VFXTs are difficult to discovery because 
their peak X-ray luminosity is below the sensitivity limits of 
all-sky X-ray instruments in orbit and therefore, most out- 
bursts are missed. Often only using pointed observations of 
more sensitive instruments, those outbursts can be detected. 
But usually, those pointed observations have a very narrow 
field-of-view and mostly are pointed towards the Galactic 
center and the Galactic bulge. Consequently the interstel- 
lar absorption is usually rather high for those systems with 
values between 10 22 cm" 2 to 10 23 cm" 2 . With such high 
column densities and the expected low surface temperature 
of the neutron stars in VFXTs, it is very difficult to im- 
possible to detect the thermal component. Even more so if 
also the non-thermal component discussed above is present 
in the quiescent spectrum as well. 

Despite that the situation described above looks very 
bleak, it might still be possible to use quiescent VFXTs 
as tests for the deep-crustal heating model. As described 
in Section 14.21 many X-ray transients are expected to be 
present in Galactic globular clusters among which many 
VFXTs (three sources are already currently known). Typi- 
cally for many clusters the column density is not very high 
and the distances are typically quite accurately known. Con- 
tinued monitoring of those clusters with sensitive X-ray in- 
struments would be crucial to catch VFXTs i n outburst so 
that later th ey can be studied in quiescence (|Heinke et alj 
l2009bl . l20ich . 

When in outburst, one specific observable property can 
be very useful to determine whether a VFXT is a primor- 
dial system or that it has accreted at significantly higher rate 
in the past: detecting millisecond X-ray pulsations. As dis- 
cussed in Section ft.ll in the context of the accreting millisec- 
ond X-ray pulsar SAX J1808. 4-3658, detecting such millisec- 
ond X-ray pulsations basically guarantees that the neutron 
star has accreted enough data to replace the crust. There- 
fore, those systems cannot be primordial (hence we expect 
spin periods >10 ms for primordial VFXTs) and the neutron 
stars systems in those systems cannot have a hybrid crust 
(unless in some way the matter spins up the neutron star but 
is not eventually accreted on the neutron star itself). The 
fact that the VFXT NGC 6440 X-2 is a millisecond X-ray 
pulsa r with a spin frequency of 4.8 ms |AItamirano et all 
l2010t ) indicates that this source is not a primordial V FXT 
and that is faint quiescent emission (|Heinke et alj|20ich can- 
not be due to non-standard heating in a hybrid neutron star 
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crust because the crust in this system should be fully re- 
placed and a different explanation (e.g., enhanced core cool- 
ing) is needed. Detailed calculations have to be performed 
whether other observable properties during outbursts can 
discriminate systems with a fully accreted crust with those 
which have a hybrid crust (e.g., type-I X-ray bursts behav- 
ior; super-burst behavior). 

Another possibility is to study those VFXTs which 
are active for a very long time (years to decades, instead 
of weeks to months), the so-called quasi-persistent sources, 
when their outbursts turn off. Severa l such quasi-persistent 
VFXTs have been identified (see, e.g jDel Santo et al.1120071 ; 
iDegenaar fe Wiinandsll2010l , albeit that those transients are 
close Sgr A* and therefore have a large column density mak- 
ing them unsuitable for studying the soft thermal compo- 
nent). It has been found that for normal transients which 
accrete so long that the neutron star crust is heated con- 
siderable out of thermal equilibrium with the core and for 
several months to years after the end of their outburst, the 
observed quiescent temperature tracks the thermal evolu- 
tion of the crust instead of tha t of the core (until thermal 
equil i brium is reached again; see Wijnands et al.ll2001 



2001 12004|; iRutledge et al.ll2002l; ICackett et al.ll200e 



2002 



2008 



201C ; iFridriksson et alj l20ld. 1201 ll; IDegenaar fe Wiinandi 
12001 



2011. 120091 : IDegenaar et alj|2009t l2011a llbl; iDi'az Trigo et al 



201 if ). If a quasi-persistent VFXT can be studied right after 



the end of its prolonged outburst, the observed crust cool- 
ing curve could give significant insight into the heating and 
thermal processes which occurred in the crust during the 
outburst. Very likely a fully accreted crust will exhibit a dif- 
ferent cooling profile than an hybrid crus10. Although this 
assumes that still a significant amount of heat is produced 
in the crust during outburst to elevate the thermal emis- 
sion above the sensitivity limits of the current generation of 
X-ray instruments. 

Obtaining crust-cooling curves for neutron stars in 
HMXBs (i.e., after the bright and extended type-II out- 
bursts observed from several neutron-star Be/X-ray tran- 
sients) will also be very interesting because it is expected 
that those curves could deviate significantly from those ob- 
tained so far for the neutron-star LMXBs. Both the effects 
of hybrid crusts in the neutron stars in HXMBs as well as 
their high magnetic field strength might produce observa- 
tional effects on the crust cooling curves. 

Primordial VFXTs will have different companions than 
systems which are only in a VFXT phase but normally they 
are normal transients or they used to be normal transients 
in the past (|King fe Wiinandsll2006l ). So finding the optical 
or IR companion star of VFXTs would be very useful to 
pin-pointed the most likely accretion rate history of those 
sources. Again globular clusters might be the best places to 
study this because of the low absorption for many of them 
compared to field targets. However, the companion star 
for both the primordial and the non-primordial VFXTs 



6 We note that for the bright transient IGR J17480-2446 in 
Terzan 5 (which might harbor a neutron star with a hybrid 
crust; see Sectionl4.H the crust cooling curve has be e n mea sured 
llDegenaar fc Wiinandsll201ll; IDegenaar et alll2011bl l2012ch and 
it did not meet the theoretical expectati ons. This could be due 
to a neutron star with a hybrid crust (see IDegenaar et al.ll2012cl 
for an in-depth discussion). 



can be very faint l|Heinke et alj|2009bl . I2OI0I ) and detecting 
them will remain a challenge, let alone obtaining spectral 
observations to confirm the type of companion star. E-ELT 
using adaptive optics might be able to detected more 
systems in the crowded globular clusters and it might be 
able to take spectra of the brightest targets. 
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